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Abstract: We study the McKean-Vlasov equation on the finite tori of length scale L 
in d-dimensions. We derive the necessary and sufficient conditions for the existence 
of a phase transition, which are based on the criteria first uncovered in [13] and [20]. 
Therein and in subsequent works, one finds indications pointing to critical transitions 
at a particular model dependent value, of the interaction parameter. We show that 
the uniform density (which may be interpreted as the liquid phase) is dynamically 
stable for 8 < and prove, abstractly, that a critical transition must occur at 9 = 9'. 
However for this system we show that under generic conditions - L large, d > 2 and 
isotropic interactions - the phase transition is in fact discontinuous and occurs at some 
9i < 9 . Finally, for H-stable, bounded interactions with discontinuous transitions 
we show that, with suitable scaling, the 9j(L) tend to a definitive non-trivial limit as 
L — > oo. 

1 Introduction 

This paper concerns the McKean-Vlasov equation - so called in [33] - which is a 
non-linear diffusion equation the classical rendition of which reads 

p t = Ap + 6L d V ■ pV(V * p). (1) 

In the above p = p{x,t), we take x G T d L - the d-dimensional torus of scale L 
- and -k denotes convolution. It is noted that the above dynamics is positivity and 
L 1 -norm preserving thus p(x, t) has a probabilistic interpretation which we relate 
to particle density. It is hereafter assumed that p integrates to one. As is well 
known, the dynamics in Eq.(l) is governed by gradient flow for the ("free energy") 
functional 

Fe(p)= [ p\ogpdx + \0L d [ V{x-y)p(x)p(y)dxdy. (2) 
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In particular all steady state solutions of Eq.(l) must be stationary points of the 
functional in Eq.(2). These densities satisfy an Euler-Lagrange equation, namely 

e -eh d [pW]{x) 



sometimes known as the Kirkwood-Monroe equation [20]. (The above follows 
from the fact, readily checked, that the dynamical equation can be recast into the 
form 



The volume factor, L d , associated with the coupling strength in Eqs.(l)-(3) may 
appear unfamiliar to some but it is in fact a principal subject of this note. 

We shall not digress with a detailed discussion of the motivations for the study 
of Eqs.(l)-(2). It is sufficient to mention the following: 

• Eq.(l) can be realized as the large iV-limit of the iV-particle Fokker-Planck 
equation under suitable rescaling of the interaction. This goes back to the original 
derivation by McKean [26] and, even today, is an active topic of mathematical 
research. A partial list of relevant papers: [5], [28], [30], [29], [3], [23]. 

• Eq.(l) can be realized as a diffusive limit of the standard Vlasov-Fokker 
Plank equation. Cf. the derivation in [24]. 

• The model of chemotaxis introduced by Keller and Segel [19] is, in fact, 
precisely the McK-V equation in slightly disguised form with a Newtonian (loga- 
rithmic) interaction; cf. [29], [15] for a derivation from particle dynamics. For our 
purposes, the Keller-Segel form of the interaction is overly singular - by no means 
a requirement dictated by biological applications. Related models with biological 
applications are described in [22], [10], [8], [7], [1], [2]. The latter two are exactly 
the McK-V equation without the diffusive term. 

• In a number of older works, beginning with [16] and [17] and including (but 
not limited to) [18] [12], [13], and [21] (cf. the article [32] for additional infor- 
mation and references) the van der Waals theory of interacting fluids in statistical 
equilibrium was elucidated as the limit of "realistic" systems under scaling of the 
interaction range. A modified version of the functional in Eq.(3), evaluated at 
its minimizer constitutes the free energy for these (limiting) theories. Finally, in 
the remarkable work [20] - predating all of the above by over two decades - the 
equation (3) for the equilibrium "distribution function" was inferred, under certain 
approximations, by direct considerations. 



p{x) 
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It should be remarked that the scaling limits achieved in the first item are not 
always in accord with those of the last. As such the volume factor is conspicuously 
absent in many modern mathematical treatments of these and related problems. 
However, on careful examination, the latter derivations contain the former in the 
static cases. Thus, for physically motivated stable interactions (which will be dis- 
cussed in Section 4) with sensible thermodynamic limits, this factor indeed belongs 
as written in Eqs.(l) - (3). For unstable interactions - which may have biological 
applications - the correct nature of the scaling has not been elucidated. However it 
appears that mathematically tractable problems in large or infinite volume emerge 
if the factor of L d is omitted. 

1.1 Mathematical Assumptions and Notations 

Since the majority of this work takes place in fixed volume, we will omit, whenever 
possible, the L-dependence in our notation for the various classes of functions etc. 
that we employ. In particular all L p -norms on T d L will be unadorned. 

The class of potentials that we consider in this work are described as follows: 
Foremost we shall assume that the V are finite range, that is 

V(x) = if \x\ > a. 

We will always take L > a and thus we may define the remaining (minimalist) 
properties as though V is a function on R d . First, we take V S L 1 and, second 
we assume that V is bounded below. The former is obviously required in order 
to make (good) sense of the uniform state. As for the latter, if V — > — oo it is 
unreasonable to suppose that this happens anywhere besides the origin. Even mild 
divergence (e.g., logarithmic in d = 2) can cause the functional to be unbounded 
below (and, in fact, just having V < a.e. in a neighborhood of the origin leads 
to unphysical behavior). Finally, on physical grounds, we shall assume that V is a 
symmetric function of its argument: V(x) = V(—x). We shall denote the class by 
Y: 

y = {V G L 1 s.t. V~e L°° and V symmetric with V(x) = for |acr| > a} (4) 

where V~ denotes the negative part of V and a < L. Additional technical assump- 
tions will be implemented as needed. 

For the analysis of the functional Tq, we shall denote by 8? the class of prob- 
ability densities on (although it is clear that 8? is much larger than necessary). 
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The uniform density, will be denoted by po : 

Po ■= L~ d . (5) 

We denote the separate pieces of Te{-) by an S and £: =: <S(-) + 

\0p^ £(■, ■). For the second, we will often have occasion to regard as this as 
the bilinear functional 

S(p a ,p b ):= V(x - y)p a (x)p b (y)dxdy 

which is (usually) defined regardless of the signs or normalization of its arguments. 
Likewise, we will have some occasion to utilize the functional S for arguments 
that, albeit non-negative, may not be normalized. For a legitimate non-negative, 
normalized p(x), the quantities S(p) and ^8L d £(p, p) are (modulo signs) vaguely 
related to the entropy and energy of the system when the equilibrium density is p; 
the two terms will indicated by these names. 

1.2 Summary and Statement of Results. 

The central purpose of this note is the study of these systems as 6 varies. Often 
enough these systems go from a quiescent (gaseous) state where no minimizers of 
Te{-) exist save for the uniform state to a state where this is no longer the minimizer 
and other minimizers are prevalent. In short, a phase transition the nature of which 
we shall partially elucidate. The results proved and their relevant location are as 
follows: 

In Section 2, the subject of phase transitions in the McK-V system will be 
discussed from the ground up. First, in § 2. 1 (which may be omitted on a prelim- 
inary reading) we establish the existence of minimizers. This allows, in § 2.2 a 
"thermodynamic" definition of the entropic and energetic content of the system as 
a function of the interaction parameter 6 which in turn will clarify the definition 
and possible nature of the (lower) transition point. In § 2.3, necessary and suf- 
ficient conditions (on V) are established for the occurrence of a phase transition. 
The candidate transition point, much discussed in other works and here denoted 
by 6 1 " is elucidated and it is shown that for 9 < 9\ the uniform density is dynam- 
ically stable. In § 2.4, a concise definition of a (lower) critical transition point is 
provided. First it is demonstrated (under the additional and presumably unneces- 
sary assumption that V € L 2 ) that if such a transition occurs, it must take place at 
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9 = OK Then it is shown that the features of a non-critical transition (where the 
above mentioned criteria fail) are dramatically different. The subsection ends with 
a principal result of this note. Namely under the majority of physically - or for 
that matter biologically - reasonable circumstances, it is a non-critical transition 
which occurs in the McK-V system. Moreover these occur at parameter value 9 T 
which is strictly smaller than 9k Finally in Section 3, the limiting behavior in large 
volume is discussed. In §3.1 it is shown that, for fixed interaction, the L — ► oo 
limit of the transition points always exist. But the limit may be trivial. In § 3.2, a 
criterion closely related to H-stability is introduced and it is shown that (with the 
scalings featured in Eqs.(l)-(3)) for stable potentials the transition points tend to a 
definitive non-trivial limit. Conversely, in § 3.3, the complimentary - catastrophic 
- cases, are investigated and it is shown that the transition values tend quickly to 
zero. 

2 Phase Transitions 
2.1 Minimizing solutions 

The starting point in our analysis is to establish, for all 9, the existence of sta- 
tionary solutions to Eq.(l) that minimize the free energy functional in Eq.(2). The 
existence of minimizers for these sorts of problems has a history: In particular [4] 
discuss the existence of minimizers for functionals of this form referring back to 
the works [13], [12]. In [6] there is an explicit construction for a related problem 
and, recently [9], established the desired result by methods not dissimilar to those 
presently employed. We shall include a proof for completeness which is succinct 
given the following: 

Lemma 2.1. Let J-q{-) be as described. Then 3Bq < oo such that for all p € 
the following holds: if\\p\\oo > Bq then there is another p$ E & for which 



Proof. We start with the observation that for any p, 



S(j>) > S(po) 



\ogL d 



and 



l -BL d £{p,p)>- l -9L d V Q 
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where — Vq is the lower bound on V{x). 

For p € and B > let B B (p) denote the set 



M B (p) = {xeT d L \p>B} 
and e 5 (p) the p-measure of B#: 

£ B = / P^X. 

We shall (rather arbitrarily) divide into the two cases of (p, B)'s for which £b{p) > 
i and e B (p) < | 

Obviously if e B (p) > 5 tnen 

S(p) > ^logB+ / plogp. 

1 Jm c b 

The second term may be estimated from below by (1 — eb) log(l — eb) + (1 — 
eb) log |B B I which can be bounded by quantities which do not depend on B. Since 
the energy term is bounded below it is clear that for some B\ < 00 if B > B\ 
Te{p) will exceed Tg{po). We turn attention to the cases eb < \ ■ 

We write p = + p r where p^ is the restriction of p to the set B b an d p r is 
the rest. Our claim is that if B is too large then 

T e {p) > Fe{p r ) 

were p r = (1 — £ B ) p r is the normalized version of p r . 

We write, S(p r ) := J Td p r log p r dx, notwithstanding the fact that p r is not 
normalized, and observe (assuming B > 1) that 

S(p) = S{ Pr ) + E B log B > S{p r ). 

Since we might as well assume that F$(p) ^ ^e(po) an d the energetic components 
of both of these quantities are bounded above and below this implies that for some 

<C 00 

> S{p) > S{ Pr ) 

regardless of the particulars of p vis-a-vis B and eb- Similarly, we have (since 
£(p, p) < £(po, po) and e b < 1/2) that £(p r ,p r ) < £{po,Po) + V =: e* < 00. 
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Now let us estimate Fo(p) — Te(p r ). First 

S{p) - S(p r ) > e B log B - -^—S(p r ) + log(l - e B ) 

1 — £b 

r i + ^i <® 

log is 



1 -£B 



where we have used the fact that log(l — eb) > — 1 f B . As for the energetics, it 
is seen that 

p) > £(p P , P P ) ~ OVqEb 

while 

£(pp,P P ) = ( 1 _ £B y £ (Pp>Pp)- 

so 

£(P,P) -£(~Pr,Pr) > [ ff B + ff ] g (Pr,Pr) > -8^* (7) 

(1 - ear 

where we have used e# < 1/2. 

The combination of Eq.(6) and Eq.(7) show that if S exceeds some (finite) Z?2 
the density p r represents an "improvement". Note that, conceivably, the improve- 
ment may take values as large as twice B2. Nevertheless, the theorem is completed 
by declaring Bq = max{i?i, 2B2, 1} and using for p* the uniform or above de- 
scribed density as appropriate. □ 

Theorem 2.2. Let J-g(p) be as described in Eq.(2) Then there exists a pg > £ 
that minimizes J~q(-\ 

Proof. Let (pj) denote a minimizing sequence for Te{-). Since, without loss 
of generality pj £ L 1 n L°°, we may place the pj in I? with | )Pj | li < -^o- 
Let poo denote a weak limit of the sequence. By standard convexity arguments, 
limj^oo S(pj) > S(poo) (where we have used (pj) to denote the sw&sequence). 
We claim that 

Km S (pj,Pj) = S (poo,Poo)- (8) 

j->oo 

This follows from some elementary Fourier analysis: Since V(x) 6 L l , \V(k)\ — > 
as k — > 00. Let uJfe = maxi^i^i where here and throughout it is 

assumed that all fc's are legitimate wave vectors for Tjf. Then 

\S{pj, Pj ) -£(p^Poo)\ < I Yl V(k)[\ Pj (k)\ 2 - \ Poo (k)\ 2 } + 2Bow ko . (9) 

fc:|fc|<|fco 
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The first term tends to zero since for each individual k, Pj(k) — * Poo{k) and the 
second term can be made as small as desired. Thus we may conclude that p^ 
actually minimizes the functional. □ 

On the basis of the above, we may define 

JCq ■= {p G & | P minimizes Fq(-)} (10) 

with the assurance that V#, 7^ 0. As an obvious corollary to Lemma 2.1, we 
have that any p € is bounded above. Conversely, we have uniform lower 
bounds (which, strictly speaking, do not play a role in later developments). 

Proposition 2.3. Let p G Then p is bounded below strictly away from zero. 

Proof. We appeal directly to the Kirkwood-Monroe equation (Eq.(3)) from which 
it is clear that pointwise upper and lower bounds on V * p are sufficient. Obviously 
V -k p < ||p||oo||^||i- Next, with more elaboration than may be necessary, let 

P a (p) = sup I p{y')dy' (11) 

yei d L J\y-y'\<°- 

where it is recalled that a denotes the range of the interaction. Then V*p > —P a Vo. 
This provides 

p(x)>e^p-eL d [P a V + \\p\\ oo \\V\\ l } >0. 

□ 

Remark It is anticipated that in physically reasonable (stable) cases, which will be 
discussed in Section 4, both terms in the square bracket appearing in the previous 
equation are of the order of L~ d . However, in catastrophic cases, it seems that 
Pa{p) will indeed achieve values of order unity independent of L for p G Mq. 



2.2 Thermodynamics for the McK-V System 

We may now separately define the energetic and entropic content of the system as 
a function of the parameter 9; these form the basis of a thermodynamic theory. 
Definition We define 

E e = inf ^9 P Q l £(p,p) (12) 
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and 



inf S(p) 



(13) 



Furthermore, defining 



inf F{p) 



(14) 



we have, to within signs and constants, the energy, entropy and free energy of the 
system at parameter value 9. It is noted that the first two do not always add up to 
the third. 

Proposition 2.4. Consider the above defined thermodynamic functions. Then 

(a) Sq is non-decreacing 

(b) Fe — ^9pQ 1 £(po, po) is non—increacing and continuous 
while 

(c) 6~ l Ee and Eg — \6p^ 1 £(po, po) are non—increacing. 

We remark that the subtractions are actually necessary: consider, e.g., the situ- 
ation - which is the general rule for "reasonable" systems - where £(po, po) > 
and JF(-) is always minimized by po for all values of that are sufficiently small. 

Proof. We shall start with the energetics. Let 0i, 02 > and let pg 1 £ ^fe 1 and 
similarly for pe 2 . Then, using pe 2 instead of pe 1 we have that 



so that (0 2 - 9i)£(p0 2 ,p0 2 ) < (0 2 - 9i)£{pe 1 ,pe 1 ) which, if 2 > 0i, certainly 
implies the first of the items in (c). However, a bit more has been shown: The en- 
ergetic content of any pe x £ is monotonically related to the energetic content 

of any pe 2 G ^e 2 ■ 

This immediately establishes monotonicity of the entropy-term. Indeed, sup- 
pose that, 02 > 0i- Then, at = 9\ using a pe 2 we have: 



Similarly, 




F dl <S(pe 2 ) + -9 lPo £(pe 2 ,pe 2 ). 



(15) 
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The energy term is less than that associated with £{pe 1 , pe 1 ) and we arrive at 

F dl < S{p d2 ) + -e x p Q - 1 £(p ei ,pe 1 ) = F 01 + S{pg 2 ) - S(p 6l ). (16) 
We again have that for any pg l € ^Kq 1 and pg 2 £ ^g 2 , with 9\ < 9 2 , 

S{p dl ) < S{pe 2 ). 

As for the claims about Fg, continuity follows from the first two displays in 
this proof. For the monotonicity, of Eg — \9p^ 1 £(po, po), we first observe that for 
any 9 and any pg € Mg 

£(pe,Pe) < £(po,Po)- 

with equality only if pg = po a.e. Indeed, assuming that pg is not a.e. equal to 
po, then S(po) < S{pg) so pg could not possibly be a minimizer if the opposite of 
the above display were to hold. Then [£(pg,pg) — £(po,po)] is non-positive and 
non-increasing so 8[£(pg,pg) — £(po, po)] is non-increasing. 

The final claim is now proved by reiteration of the previous procedures with 
the subtraction in place: 

F 02 -^6 2 po 1 £(po,po) < S(p 6l ) + ^6 2 pQ 1 £(pg 1 ,pg 1 ) - ^0 2 Po 1 £ (po,Po) 

= F 9i - 2 iPo l£ (po,Po) + \iP 2 - 9i)p 1 [£(pe 1 ,pe 1 ) - £{po,Po)} 

where we have assumed 6 2 > 0\. By the non-positivity of the quantity in the 
square brackets, the stated monotonicity is established. □ 

With the above monotononicities in hand, the objects Sg and Eg can now be 
considered well defined functions for all 9 which are continuous for a.e. 9. How- 
ever at points of discontinuity, it may be more useful to focus on the range of the 
function rather than its value at the point. In particular Sg is continuous iff Eg is 
continuous while at points of discontinuity, the density that minimizes S is the one 
that maximizes £ and vice versa. 

2.3 Phase transitions in the McK-V systems (1): 
The point of linear stability 

We start this subsection with some preliminary results - most of which have ap- 
peared elsewhere in the literature (albeit by different methods) - concerning the 
single phase regime: The regime where po is the unique minimizer of Tq{-). 
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Proposition 2.5. Let V E "¥ be bounded i.e. |V| < V max < oo. Then for 9L d 
sufficiently small — less than [Vmax]" 1 - the functional Tq(p) is convex. 

Proof. The functional F${p) is finite on the set 

£ = {p E & | plogp E L 1 }; 

and there is no ambiguity to set F$(p) = +°° for p E &\J2. 

Then, it suffices to show that for any p\,p2 E J2 the function s i— > Feips) 
where p s = /32s + pi(l — s) is convex. It is straightforward to verify that Fe{Ps) 
is twice differentiable in s E (0, 1). Then, we compute 

,2 



(P\ 2 HPs)= [ ^dx + OL 
\dsJ Jjd p s 



T i xT i 



V(x — y) r](x)rj(y) dx dy, 



where r) = p% — p\. By Jensen's inequality we have 



V-dx 



If Ps 



V 



jd \p s 



— p s dx> 



p s dx 



\rj\ dx 



On the other hand, since \V(x — y)\ < V max we have 



/ 



V(x — y) T](x)rj(y) dx dy 



This implies the inequality 

(-^) 2 Hps)>(i-eL d v w 

H6L d V max < 1. 

This immediately implies: 



\t]\ dx 



77 1 dx ) > 



□ 



Corollary 2.6. Under the conditions of Proposition 2.5, po = 1/L d is the unique 
minimizer of Tq (p). 

From the proof of the above Proposition we also obtain a corollary for poten- 
tials V which are of positive type: 

Definition A potential V is said to be of positive type if for any function h E L 1 , 



T i xT i 



V(x — y)h{x)h{y)dxdy > 0, 
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which is equivalent to the condition that Vfc 

V(k) > 0. 

We let "V + C "V denote the set of interactions that are of positive type and, for 
future reference, the complimentary set by Y^: 

y N = y\y+. (17) 

Corollary 2.7. Let V € Y + . Then for all 9, the unique minimizer ofj-g(-) is the 
uniform density pq. 

Proof. For p = po(l + rj) in we consider /^(s) := ^(p^l + sr))). Calculating 
fq(s) as in the proof of Proposition 2.5, the entropy term is still positive while the 
energy term yields 

p / V(x - y)r](x)r)(y)dxdy > 0. 

Thus f v (s) is always convex and, for all rj always minimized at s = 0. 

Note that since all convexities are strict, any p € 8? that is not a.e. equal to po 
admits, for all 9, 

Fg(po) < Fe(p)- 

□ 

Next we show that V G Yn is also sufficient for the existence of a non-trivial 
phase. The starting point is an elementary result which, strictly speaking is a corol- 
lary to Proposition 2.4. 

Proposition 2.8. Let V G and suppose that at some 0^ < oo there is a pg d 
which is not a.e. equal to pq such that 

Fe d (pe d ) < Fe d (po)- 
Then then for all 9 > 9^, Po is not the minimizer ofTg(-). 
Proof. Indeed, since pg d is not a constant it must be the case that 

S( P6d ) > S(po) (18) 
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thence 

£(pe d ,Pe d ) < £(po,po). (19) 
Thus for 9 > 84, it is seen that 

F e < Fe d (pe d ) + ^PoHO - 9 d )£{ P0d ,pe d ) 

< Fo d (pe d ) + gPo 1 ( 9 ~ Od)£(po,po) < Fe(po) 
which is the stated result. □ 

Thus beginning at 9 = there is a non-trivial region or phase characterized by 
the property that po is the unique minimizer for Te(-) and this phase terminates at 
some value of 9 - which is possibly infinite. Assuming this value is finite, we may 
refer to it as the lower transition point and above this point, there are non-trivial 
minimizers of Tq and non-trivial solutions to Eqs.(3) and (1). (We shall refrain 
from naming this point till the possible nature of the transition at this point has 
been clarified.) It should by noted, by a variant of the above argument, that at the 
lower transition point po is actually still a minimizer of the functional in Eq.(2). 

We introduce some notation: 
Definition For V G "V^s, let fe" denote a minimizing wave vector for V(k): 

t>(fe") < V(k) Vfe. 

Note that V(0) < by assumption. We define #" = 9\V) via 

0" := \V{^)\~ 1 . 

We are finally ready for the following: 

Proposition 2.9. [13]; see also [4] Let V G Y N . If 9 > 6>" then 3p G & , p + p 
which minimizes J-g{-). In particular, for 9 > 9^, po is no longer a minimizer of 
Tq. Thus V G is the necessary and sufficient condition for the existence of a 
non-trivial phase. 

Proof. For 9 > 9^ we may use as a trial minimizing function 

P = P o{l+erp) 
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where rf 1 is a plane wave at wave number k^ and is itself of order unity while e 
is to be regarded as a small parameter. Since all quantities are bounded, we may 
expand: 

po(l + er/)logpo(l + er/) = 

p (l + erf) log/*) + p (l + erf){eif - ^[eif] 2 ) + o(e 2 ). 

Since 77" integrates to zero, 

S(p) = S(p ) + ^e 2 Po J \rf\ 2 dx + o(e 2 ). (20) 

Meanwhile 

l -9p G 1 Z(p,p) = hp^£{ Po , Po ) + ^ 2 0p o / V(x - y)r,\x)rf{y)dxdx 



By definition of i 



^p -^(po,p ) + \e 2 p o V(0)\\^\\ 2 ^ + (9- 



-^ 2 lMll2 = ^ 2 ^)l^lli 



so that 



^(p) = ^(Po) - \e 2 Po\W\\l W^)\] (9 - 0») + o(e s 



which is strictly less than Tq (p) for e sufficiently small. (Here it is noted that the 
quantity polMlH * s itself of order unity.) By Proposition 2.8 the above is sufficient 
to establish the statement of this proposition. □ 

Corollary 2.10. For V G "Vn, 6'(V) is the supremum of the set of quadratically 
stable parameter values for Tg(po). Furthermore, 9^ marks the boundary for the 
linear stability ofEq. (1) with solution po I.e. for 9 < 9\ po is linearly stable while 
for 9 > 9^ it is not. 

Proof. The first statement follows, in essence, from the above display. As for the 
dynamics, the linearized version of Eq.(l) reads, for rj = (p — po)Po 1 

^=V 2 (ri + 6V*Tj). (21) 

The linear operator V 2 [l + QV*](-) has, by the definition of 6$, a strictly negative 
spectrum if and only if 9 < 9K The second statement of this corollary therefore 
follows from the definition of linear order stability. □ 



14 



While the above is ostensibly vacuous, in these cases, it turns out that po actu- 
ally has a non-trivial basin of stability for 9 < 9*. 

Theorem 2.11. Under the regularity assumption 

G = Y,\V{k)\\k\<™ 

k 

there is a non-trivial basin of attraction for pq which contains all Borel measures 
that are sufficiently close to po in the total variation distance. In particular, at 
positive times, any such perturbing measure regularizes and, for any particular 
Sobolev norm, the density converges to po exponentially fast in this norm. The 
stated results hold uniformly in L. 

Remark The regularity assumption on V is for convenience; presumably a stronger 
result is available. In particular, it is not hard to see that with greater regularity of 
the perturbing density, regularity assumptions on the interaction potential can be 
relaxed. Moreover, with greater regularity assumptions on V, even more singular 
objects than Borel measures are contained in the basin of attraction. 

Proof. We write p = po(l + 77) with rj measure valued. Let fjk(t) denote the dy- 
namically evolving k th Fourier mode. We shall assume that in the initial state, each 
mode is small - which is certainly implied by the smallness of the total variation 
distance. In particular, we will assume that at t = each ?% (0) is bounded by an 
£0 which satisfies the condition that for all k, 

2\k\9Ge < X(k) (22) 

where X(k) = k 2 (l — 9V(k)) is the decay rate for the the k th mode in the linear 
approximation. It is emphasized that A(A;) > ck 2 with c > if 9 < 9\ 
The fjk (t) satisfy the formal equation 

= -X(k)Vk(t) + 0k-J2 k'V(k')f) kl (t)f) k ,- k {t) (23) 

k' 

where, it is reemphasized, all factors of volume have canceled out. It is noted that 
Eq.(23) may certainly be used to formulate dynamics via an iterative scheme - 
provided that control is maintained under reiteration. Thus we may consider the 
sequence (77! (i) | I = 0, 1, 2, . . . ) where fj^it) = %(0) and ff^ +l \t) is defined 
as the solution of Eq.(23) with f] k (t) the argument of the non-linear kernel. 
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The form in which we will use equation (23) is with moduli; we have 

= -X(k)\fj k (t)\ + ^0[^k.^k'V(k')f, kl (t)f, k ^ k (t)+c.c. 

The first claim is that for eo satisfying the condition in Eq.(22) then for all k and t 

and t, 

\v { k \t)\<e . 
Indeed, this is certainly true for so, inductively, 

^^<-X(kU e \t)\ + 6elG\k\. (24) 
First, let us consider modes that satisfy 

h»(0)l > LLj*. (25) 

Such modes will decrease in magnitude - at least till |r/fc(i)| reaches the right side 
of the inequality in Eq.(25) whereupon they may "stick". But by assumption, these 
modes started out smaller than eo- On the other hand, modes with initial conditions 
that satisfy the opposite inequality of Eq.(25) may actually grow till the inequality 
saturates but this does not get them past eo since for all k, 

by hypothesis. (The factor of two does not yet come into play.) 

Contraction of the sequence follows an identical argument which does employ 
the factor of 2. We define A e k (t) = \fj^ +l ^ - fj^\ and = sup M Af (t). It is 
found that Vfc, t, 

Ai(() <^J^ <(1 _, K -. (27) 

for some S > by Eq.(22). Thus Eq.(23) indeed defines our dynamics and we may 
perform manipulations on its basis without further discussion. Our next task will 
be to get the r] k uniformly decaying. 

By repeating the steps of Eqs.(24) - (26) it is clear that for any eo > eo, there 
is a time to such that for all t > to, 

, . . G6e?,\k\ 

\m(t)\ < -^J. (28) 
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Incidentally, we have now placed 77 in some reasonable Sobolev space - but this 
is not yet relevant. For the moment, the pertinent observation is that there is a 
maximum sized mode which is to be found at a finite value of k (which may, of 
course, change from time to time). Thus, for each t > to, let /3q denote the modulus 
of the maximum mode and k denote the wave vector that maximizes. Then, for all 
t > to we have 

^ < -A(fc)A) + G6\k\p e < —XrrmPo (29) 

where A m i n is the minimum of X(k) (which is positive for 9 < 6*"). We conclude 
that all the r/fc(i) tends to zero exponentially fast with rate at least as large as ^A m i n . 

We use a small variant of this argument to show that for any n, the maximum 
of |A;| n |77fc(t)| (exists and) decays with a rate at least as large as ^A m ; n . Focusing on 
n > 1 let us assume that at the n — 1 st stage of the argument, we have a t n -\ such 
that for alH > t n -\, 

/T ^) < (30) 

here (3 k n \t) := |^| ri |?7fe (^) I an d the quantity 6 n is specified as follows: Multiplying 
both sides by \k\, since A > ck 2 this puts a uniform bound on (3 k n \t) which is 
stipulated to be less than one. 
We now write 

< -A(fc)/t ] + 2^\k\6^pP\ m/ . k \\V(k')k'\ + 

k' 

2 n - 1 \k\eY / ^-k\vk>\\v(k')k>\ 
k' 

where we have used \k\ n = \k' + k - k'\ n < 2 n - 1 [|/c| n + \k - k'\ n ]. We now 
wait till a time t' n when each \rjt(k) \ is less then some e n which is small and to be 
specified. Summing, the estimates, 



< -\{k)p[ n] + 2 n 9\k\e n (31) 



dt 

and we now wait till a time t n > t' n so that, similar to the previous portion of the 
argument, 

Pt\t)<°^^n (32) 
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for any S n > s n . Obviously this 5 n will be tailored to satisfy the requirements to 
propagate the next iterate of the argument - and so we stipulate. But in addition we 
require that 2 n 5 n < Eq 

The proof is completed by noting that Eq.(32) allows us to conclude that the 

In] 

supremum of pi is to be found at a finite k and the rest of the argument proceeds 
as described in the vicinity of Eq.(29). The desired result has been proved. □ 

2.4 Phase transitions in the V-McK systems (2): 
First order transitions. 

In order to investigate the possibility of continuous/discontinuous transitions in this 
model, an appropriate definition must be provided. 

Definition Consider the V-McK functional Tq with V £ ^n- We define 9 C to be 
a (lower) critical point if the following criteria are satisfied: 

• For 9 < 9 C , po is the unique minimizer of Tg^). 

• For 8 > 9 C , 3p0 ^ po which minimizes Tg{-) 

(which necessarily implies that po is no longer a minimizer of Tq{-)). 

• If (pq I 9 > 9 C ) is any family of such minimizers then 



Remark We have called this a lower critical transition since, conceivably there 
could be later (in 9) transitions of this type with non-trivial solutions "bifurcating" 
from preexisting non-trivial solutions. This would be difficult to detect - analyt- 
ically or numerically - since the non-trivial solutions are anyway evolving with 
9. Such a phenomenon would, presumably, have to be tied to non-analyticity in 
£ (■) or S(-) notwithstanding their continuity. By contrast (c.f. Proposition 2.13 
below) for the other possible type of transition, these objects are generically dis- 
continuous. In any case, the foremost possible phase transition in these systems is 
the lower one and will be the focus of all our attention. 

Any (lower) phase transition not satisfying the above three items will be called 
a discontinuous transition and we will denote will denote such a transition point 
by 9 T . As we shall see later, in Proposition 2.13, for a discontinuous transition the 
second item will hold while in the first item, we must replace 9 < 9 C with 9 < 6 r . 
But most pertinently, the third item fails in its entirety. Thus, at such a transition 
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point, a new minimizing solution of Eq.(l) appears which, for 9 = 6 T , is degenerate 
(in the sense of minimizing with po but is markedly separated from po- 

Our first result characterizes the critical transitions: 

Proposition 2.12. Let V £ YnHL 2 and suppose a (lower) critical phase transition 
as described above occurs in the V-McK system at some 9 C . Then, necessarily, 



Proof. The trivial cases 6$ = or 6> s = oo are easily dispensed with. Assuming 
otherwise for 9$, it is obvious that a (lower) critical point 9 C could not exceed 9* 
since non-trivial minimizers already exist at any 9 > 9$. We shall therefore work 
with 9 < 9^ and write 9 = 9$ - 6 where 5 > 0. 

As a preliminary, it should be noted that while the third item in the definition 
of the 9 C necessarily reflects the natural Z^-norm, it will be more convenient to 
work with L 2 and L°°. We will show that, as far as pe — po is concerned, these are 
controlled by the Z^-norm. First, for expositional ease, let us define 



Starting with L 2 , recall from the "obvious corollary" to 2.1 that since pq is a min- 
imizer of Tq{-) it is bounded uniformly (enough) in 9 and thence rjg is similarly 
bounded - say by uj. Then 



For the moment, we can only employ the outer inequality but at least we now have 
that ||% lb is "small". Next we use the fact that pq satisfies the Kirkwood-Monroe 
equation, Eq.(3). As is not hard to see, in the language of t]q this reads 



9 C = 9k 




(34) 



e 



[9V*r,g](x) 



1 + rje(x) 



(35) 



Now, for a.e. x 



\{V* Ve }(x)\= [ V{x-y)i le {y)dy <\\V 



2\\f]e\\2 



(36) 



thence, if r)g(x) > 0, 



Ve(x) < (e 2 *!!^ 2 ^!! 2 - 1) 



(37) 
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while if f]e(x) < 0, 

Ve(x) > ( e - 2e \\ v \\2\\ve\\2 _ X ) (38) 

Thus, for \ \rjeW2 sufficiently small (which we know happens as 9 J, 9 C from Eq.(34)) 
there is a K - which is uniform in 6 near 9 C and of order unity - such that 1 1 ^76> 1 1 00 < 
K\ \no\ b- We may now exploit the middle inequality in Eq.(34) and declare that in 
the vicinity of the purported 9 C all norms of any rjg are comparably small. 

Now, suppose that 9 > 9 C . We repeat the calculations performed in Proposition 
2.9 with the result that 

Fe(pe)=S(po) + ^9L d £(p , Po ) + ^p [\\ Ve \\ 2 2 + 9£(m,Ve)] (39) 

The term in the square brackets is strictly positive and at least of the order \\rig\\\ 
if 9 « 9 C = 9^ — 5 with 5 > 0. Evidently, as indicated, the only possibility for a 
continuous transition is at Ok □ 

The alternative to a critical transition is a discontinuous transition which is also 
called a first order transition. For such transitions, the following holds: 

Proposition 2.13. If V £ "Vn and the criteria in the preceding definition of a 
critical point fails then there is a transition at some 9 T which is characterized by 
the following: 

3pe T 7^ po such that 

• Fe T {p8 T ) = Fe T (po) = Fq t 

• £(po T ,Po T ) < £{po,Po) 

• S{pe T ) > S(p ) 

( and thus both Eg and Sg are discontinuous at 9 = 9 T ). 

Since two distinctive minimizers exist at the same value of 9, such a point may 
also be described as a point of phase coexistence. 

Proof. At 9 > 9 T we have for r]g = (pg — po)Po 1 

limsup ||??6»||i / 0. (40) 

Since, in these matters, all norms are more or less equivalent, we will take the 
above statement in I? and extract a weakly convergent sequence which we will still 
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denote by r\g. Let us first rule out the possibility that rjg — 1 0. Indeed, supposing 
this to be the case, we would certainly have 

lim £(rje,r)g) = 0- 

v — >UJ 

e.g., as discussed in the proof of Theorem 2.2. However, we have that all along 
the subsequence, \\r]e\\2 > h T for some h T > and, moreover, for some b < oo, 
1 1 770 1 1 oo < b. Thence, by the convexity properties of the 5-term we have that for s 
small, 

S(po) >S(p ) + ^s 2 \\ V e\\ 2 2 + o(s 2 ). 

This indicates that 

lim sup Fg > Fg T 

in violation of the stated continuity result. 

Thus, in our sequence r]g converges to a non-trivial limit which we denote 
(optimistically) by r]g T . On the energetic side, we still have 

lim £(rjg,rj e ) = £(r]g T ,r]g T ) 

6 — >uj 

and, again, by convexity properties, S{p${l + r]g T )) does not exceed any limit of 
S(pg) as 9 I 6 T . Evidently this ?]g T provides a genuine minimizer for Tg T {-) which 
we now denote by pg T . 

By hypothesis (of a lower transition) the uniform solution is a minimizer of Tg 
up to 9 = 9 7 and thus by continuity is also a minimizer at 9 T : Tg T (po) = Tg T {pg T ) 
(see Proposition 2.4). Moreover, we reiterate, <S(po) < S(pe T ) necessarily imply- 
ing £(pg T ,pg T ) < £(po, po). All of the stated results have now been proven. 

□ 

The two preceding results - concerning (i) the purported critical behavior at 
9 = 9$ and (ii) the characteristics of systems with purported non-critical lower 
transitions - allow for the following: 

Theorem 2.14. Consider, in dimension d > 2 a fixed V € 'Vn which is isotropic. 
Then, if the volume is sufficiently large there is never a (lower) critical transition. 
In particular under the above stated conditions there is a discontinuous transition 
at some 9 T satisfying 9 T < 9$ where there is phase coexistence and various other 
properties all of which has been described in the context of Proposition 2.13. 
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Proof. We will consider disturbances of the form 

p = p (l + er)) 

with ?](x) a function (with L°°-norm) of order unity and e a small (pure) number 
of order unity. Then 

S(p) = S(po) + / P v[\eW - \e*r?]dx + o(e 3 ) (41) 
J T d 2 b 

where it is slightly important to observe that o(e 3 ) is independent of the volume. 
Of course the above expansion also contained a linear odd term which vanishes 
due to symmetry. Similarly, we have £(p, p) — £(po, po) = e 2 £(por], porj). 

We set = 6^ where, as we recall, the minimizing wave vector satisfies 
— V(k^)6^ = 1. Now, we invoke the assumption that V{x) depends only on 
\x\, - so that V(k) depends only on \k\. Then, under the auspices of contin- 
uous wave numbers ("the infinite volume limit") we could find k\ and ^2 with 
\k^\ = I k\ I = I ki I necessaiily satisfying 

V(0) = V(h) = V(h) (42) 

such that 

k i + k 1 + k 2 = 0. (43) 

Thus, in finite volume, we can find approximating k\ w k\ and k 2 ~ ki with, e.g., 
\k\ — k\ \ = OiL" 1 ) that are appropriate to such that Eq.(43) is true and Eq.(42) 
is approximately true. We now use 

77 = rf + 7/1 + r}2 

with 7/1 and 7/2 plane waves at wavenumbers k\ and k 2 respectively. We have, e.g., 

Vihppol^ih)^ + p \\ m \\ 2 2 < a(L) (44) 

with a — > as L — > 00. (We reiterate that each term in the above display is sep- 
arately of order unity.) Thus, we may declare that, essentially, up through second 
order !Fgt(po(l + erj)) equals Fg^po). But now, since fe" + k\ + k 2 = 0, then 
unlike a plane wave which, even cubed, would integrate to zero, it is in general the 
case that 

/ (r] 3 )dx f 0. (45) 
hi 
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Moreover by adjusting the phases of the constituents, the corresponding term in 
the expansion of T e i (p) can be made to be negative. 

It is thus evident that for all L large enough, po is not the minimizer for T e % (•) 
and by the continuity result contained in Proposition 2.8, po is n °t the minimizer 
for a range of 9 which lies strictly below dk Thus the transition takes place at some 
9 T < 6$ and is (therefore) not continuous. 

□ 

We conclude this section with an (abbreviated) spectrum of remarks. 

Remark For the vast majority of physically motivated single component systems, 
the above theorem precludes, in the general context, the possibility of continuous 
transitions. (Cf. the third remark in this sequence for additional discussion.) This 
is in apparent contradiction with a number of results for these system - some of 
which receive additional discussion in the subsequent remark - the most pertinent 
of which have been the subject of [6] and, recently, discussed in [4] . In these 
works, a continuous transition was indeed found at the analog of 9k The important 
distinction distinction between the present work and [6], [4] is in the nature of the 
entropy functional that was employed. Indeed, therein the prototypical entropy 
functional was of the form 

MP)= f[p\ogp+(l-p)\og{l-p)]dx. (46) 

Thus, in the expansion which uses p = po(l + r/), all the odd terms in 77 vanish 
identically; from this perspective, Ao(p) is simply the symmetrized version of 
S(p). Of course this preempts the term(s) driving the conclusion of Theorem 2.14 
and thus allows for a continuous transition at 9 = 8\ 

However, Ao(p) is not a natural entropy form for a one-component system 
and, as argued in [4], is in fact an effective entropy term for a two component 
Ising-type system. The first principles version of these sorts of Ising systems is 
currently a subject of intensive investigation; e.g. the works [6] and [11] and some 
work in progress by R. Esposito and R. Maria in conjunction with the authors. The 
present set of models under consideration appear to undergo a transition that is, 
at least sometimes, "weakly first-order" at some 9 T < 9k However, it may well 
be the case that the consideration of more general interactions leads, in the two- 
component cases, to generic circumstances where there are continuous transitions. 
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Remark In a variety of contexts, e.g. [14], [31] [25], various workers have 
claimed that non-trivial solutions to Eq.(3) grow - continuously or otherwise - 
only after 9 > 9\ Since this attitude seems prevalent, perhaps some comments 
are in order. As is typical - by definition - for discontinuous transitions the new 
minimizing solutions or states are not continuously connected to the old. This and, 
especially, Theorem 2.11 accounts for some of the difficulties attempting to gener- 
ate stable solutions dynamically as described in [25]. It may, perhaps, prove useful 
to attempt to generate the stable solutions by a nucleation technique (e.g., based on 
the solution for small systems) at parameter values below 0" where, perhaps, fewer 
interfering solutions exist. Indeed, the basin of attraction provided by Theorem 
2.11 may itself, for L ^ 1, be arguably small from a certain perspective. 

The results in [14] and [31] both rely on standard fixed point/bifurcation anal- 
yses. In [14], it was simply assumed that the non-trivial solutions were periodic 
with period 0. Of course, as was noted in [14] the Kirkwood-Monroe equation is 
"closed" under periodic functions with any period. (By this it is meant that if we 
write Eq.(3) in the form p = H(p) then, if a is periodic so is H(<r) with the same 
period.) Thus, using this equation as the basis for a fixed point argument (with 
the help of the Krasnoselskii fixed point theorem) one is liable to manufacture a 
solution of sorts. Moreover, this scheme indeed requires 9 > 9^ for the solution of 
period 27r[fc"]~ 1 to be non-trivial. However it is also clear that these solutions have 
no stability under the dynamics of Eq.(l) and, even, the discrete time dynamics 
which produced these solutions in the first place. In particular it is almost certain 
that these solutions do not minimize the free energy functional. (Although, no 
doubt, they have a lower free energy than the uniform solution.) Indeed while it is 
not impossible that the stable solutions appearing at 9 = 9 T are periodic with some 
period, as is typical in non-linear phenomena, there is no reason for the period to 
exactly match that of the unstable mode which appears at 9 = 9$. 

Finally, we wish to comment on the careful analysis in [31]. Here, standard 
results in the theory of bifurcations were brought to bear under the explicitly stated 
assumption that the relevant hypotheses for the theorem actually apply. In this 
context, the most important of these ingredients is that the kernel and co-kernel of 
the linear operator are one-dimensional. Under the required symmetry V{x — y) = 
V(y — x) - without which the model does not make sense as a description of 
identical interacting constituents - this condition is obviously violated. And it 
may or may not be a "technical" violation, cf. the next remark. Notwithstanding, 
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even if the conditions for the bifurcation results are satisfied, provided that V(k) 
is continuous, there are always modes near fe" which are nearly unstable at = OK 
Thus any basin of stability and domain of validity will be vanishingly small with 
increasing volume. 

Remark It is remarked that full isotropy of V(x) and/or d > 1 is not strictly re- 
quired. The condition used in the preceding proof was the existence of three wave 
vectors adding up to zero each of which (nearly) minimize V(k). Obviously this 
can be achieved in d > 1 if V(x) has an appropriate 3-fold symmetry. Moreover, 
a detailed analysis will yield alternative sufficient conditions: (I) If V"(0) (assumed 
positive) is not too large. (II) if V(20) is negative and, in magnitude, an appre- 
ciable fraction of V(0); etc., etc. However, full isotropy is not an unreasonable 
assumption for fluid systems - as well as other applications - and, in fact, d > 1 is 
required for actual statistical mechanics systems with short-range interactions to 
exhibit changes of state. Thus we are content with the present result and will not 
pursue these alternative specialized circumstances. 

Remark In the language of equilibrium statistical mechanics, 6 T is, of course, 
classified as a point of first order transition while the point 9^ is not recognized. 
From the perspective of dynamical systems, 6^ is a subcritical pitchfork bifurcation. 
It may be presumed that solutions of the type which minimize at and above 9 T 
are present even before 6 T . The point at which they first appear - temporarily 
denoted by #r - would then represent a saddle node bifurcation while, from this 
perspective, the point 9 T is not recognized. 

3 The large volume limit 

In this final section, we shall investigate the behavior of our systems - with fixed 
V(x) - as L tends to infinity. The upshot, roughly speaking, is that for interac- 
tion potentials which are appropriate for physical problems the energy/temperature 
scaling is viable and not so otherwise. Since the L-dependence of these problems 
will now be our focus, all relevant quantities will adorned with superscript [L\. 

3.1 The limit of the transition points 

If the interaction violates the conditions of Theorem 2.14 and has a (sequence of) 
continuous transitions then, by Proposition 2.12, these all take place at the relevant 
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9^ which is only weakly dependent on system size. Thus the interesting questions 
concern the discontinuous transitions. Notwithstanding, the forthcoming makes no 
explicit use of the discontinuity other than the convenience of label. 

Theorem 3.1. For fixed V € Vn consider the system on with discontinuous 
transition at 6j. Then these transition points tend to a definitive limit. 

Proof. We shall start with the statement that for any L and any integer n, 

d\ L] > e[ nL] . (47) 

To see this, we patch together n d copies of a non-trivial minimizer of T d L at 9^ 
to cover T d L (which is facilitated by the fact that, anyway, these solutions are 
periodic). First, letting 

ix (48) 



/ V(x)d2 

jR d 



it is noted that for any L a , 



4 La] (Po a] ) = -^gL d a + hv. (49) 

Now let L denote any scale with transition temperature 9^ and let denote the 
non-trivial minimizer for at this value of the parameter. Let pi™ 1 ' denote the 
periodic extension of this function to rescaled by a factor of n~ d so that it is 
properly normalized. It is seen that 

S nL {p [nL] )= I' p [nL] logp [nL ]dx 
JT d T 

= - log n d + n d x 4j / pl L] log pf ] dx 
n d Jjd 

= -logn d + S^(p[ L] ). 

Making use of the underlying periodic structure, we have that for fixed L- 
periodic g(y), the integral J T<i V{x — y)g{y)dy is equal to the periodic extension 

nL 

of the corresponding integral on T^. Thus, the energetics will come out the same. 
In particular, if we define 

N(x)= [ p [ : L] {y)V{x-y)[nL] d dy (50) 
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then N is the periodic extension of the function N(x) given by 

N(x) = [ p[ L] (y)V(x-y)L d dy (51) 

(Note that in Eqs.(50) - (51) the factors of (nL) d and L d have been brought inside 
so that the integrands are both ostensibly of order unity and therefore the "same" 
function.) 

Thus, again, 

\e[nL] d E^{pt L \p [ : L] ) = \en d x ^L^ L \pf\p^). (52) 
Altogether, we find 

^ L \p^) = -logn d + F [ e L] (p[ L] ) (53) 
while, from Eq.(49) with L a = nL 

4 nL] (Po L ) = -logn d + F [ e L] (p [ L] ). (54) 

From the above two equations, we may conclude that #[ L ' > 6^ L \ 

Now consider L, K with K 3> L when K is not an integer multiple of L. We 
will use almost exactly the above argument except that we will acquire an error due 
to "boundary terms". Let us find n such that 

nL < K < (n + l)L (55) 

We shall treat T^- like the hypercube [0, K] d which is divided into n d hypercubes 
of scale L which occupy [0, nL] d . For future reference, we refer to boxes that share 
a face with the region T^- \ [0, nL] d as boundary boxes. It is noted that there are 
B(n, d) = n d — (n — 2) d such boxes. In the region [0, nL] d , we define, similar to 
before, the density p^ which is the rescaled periodic extension of p^ , the non- 
trivial density which minimizes the free energy at o[ L ^ on T^. In the complimentary 
region, we set p\ to zero. 

The entropic calculation proceeds exactly as before with the same result namely 
— \ogn d + S^(p^). However, for the energy integrals, we cannot simply use 
Eq.(52) because, e.g., if both x and y are in boundary cubes (on opposite sides) the 
formula may be in error because V(x — y) no longer "reaches around". However 
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for present purposes, it is sufficient to use the results of Eq.(52) and subtract the 
maximum possible gain from these cubes - which would be —Vq. The result is 

rt(pFi,^)<(^)V^a L] ,^ ] )]+ 

+ V -[nL] d [ pf\x)pf\y)dxdy] 

JH(n,<i) 2 

where H(n, d) is standing notation for the above described region of boundary 
cubes. Note that the p^'s are normalized to n~ d in each such cube and the order 
of B(n, d) is n d ~ 1 . As a result, 

K d £{pf\pf ] ) < [L d £(p[ L \p[ L] )](l + 0(K/L)) (56) 

Consequently, for any e > 0, 

ef ] + e > e[ K] (57) 

for all K sufficiently large which implies the desired result. □ 

3.2 Stable Behavior 

Since o[ L ^ tends to a definitive limit which (for V £ Yn) is not infinite, it is im- 
portant to establish the criterion for when this limit is not zero. As it turns out, the 
correct condition is closely related to thermodynamic stability. 

For two body interactions, the condition of H-stability (see [27] p. 34) is as 
follows: 3b > — oo such that for any N points, x\, . . . x^ in R d , 

^V{xi-Xj)>-bN. (58) 

This condition is necessary and sufficient for the existence of thermodynamics - 
although the early proofs usually assume continuity properties of the interaction. 
Provided that V is bounded and continuous, H-stability is equivalent to the condi- 
tion that for all probability measures described by a density p(x), 

[ V(x - y)p(x)p(y)dxdy > (59) 

JR d xK d 

(c.f. [27]) as is easily seen by utilizing sums of point masses to approximate prob- 
ability measures. This will be our working hypotheses for the benefit of the next 
result along with the technical assumption that V is bounded: 
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Definition An interaction V € "¥ is said to satisfy condition-K if |V(x)| < 
V max < oo for all iER^ and if for all L sufficiently large, the inequality 

f V(x - y)p(x)p(y)dxdy > (60) 
JT d L xf d L 

holds for all p £ 

The principal result of this section is as follows: 

Theorem 3.2. Let V £ "f^ denote an interaction that satisfies condition-K and 
has (for all L sufficiently large) discontinuous lower transitions at 9^ on T^. Then 
the 9^ tends to a limit that is strictly positive. 

Proof. For connivence in the up and coming we shall streamline notation - e.g., 
revert to the omission of all L's in the superscripts, etc. We start by assuming 
9 > 9j and write, for this value of 9, a non-trivial minimizer and its deviation 

p = p (l + ??)• 

Further, we define the positive and negative parts of rj as rj + and, rf~ respectively 
and, finally, 

h = HamIIi- 

The aim is to show that if 9 is small than, regardless of L, h must be zero. 
The first step will be an estimate on the free energetics. We have 

< Fe(po) ~ MP) = S(po) - S(p) + l -9L d [£{p GlPQ ) - £(p,p)}. (61) 

As has been stated before, we have ^9L d £(po, po) = \9v while 

^9L d £(p, p) = X -9v + l -9L d £{p Qr] , m ). (62) 

Let us decompose: 

hL d £(p ri,p r]) = ^9L d [£ (p r/ + , p r] + ) + £(poV~, PoV~)~ 
2£(p r] + ,poV~)}- 
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The first two terms are positive by the hypothesis that V satisfies the condition-K 
thus 

1 1 f 

-9L d £(p, p)>-6v-9 V(x - y)r]~(y)poV + {x)dxdy 

1 1 JT d T xT* 



> -0v - -9\\V\\ih 

where we have used | j 77 [ 1 00 < 1 and ||pof? + ||i = \h. 
Putting these together we have 

S(p)-S( Po ) < ^9L d [£(p ,p )-£(p,p)} < lemh. (63) 

Incidentally we may use the lower bound (See, [33] p. 271) S(p) — S(po) > \h 2 
to learn that the assumption that 6 is "small" necessarily implies that h is small but 
the particulars of this bound does not play a major role. 

Now, let us write the mean-field equation, Eq.(3), in a form useful for the 
present purposes: 

log p + 0L d [ V(x-y)p(y)dy = C KM (64) 

with Crm a constant that we are now prepared to "evaluate" 

C KM = S(p) + 9L d £(p,p). (65) 
Expressing Eqs.(64) - (65) for the benefit of rj we have 

log p + log(l + V ) + ev + e[ V(x - y)i](y)dy = S(p) + 0L d £(p, p) (66) 

so 



log(l + 7]) + 9 / V(x - y)r](y)dy =: -k 
JT d r 



Fe(fi) ~ Mpo) + \9L d [£{p,p) - £(p , Po )] 

where it is noted that the sign of k is pertinent. Indeed by the display just prior to 
Eq.(63) we have 

< k < -9\\V\\ih (67) 
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(Thus, in addition, k is small). Note that all the above holds formally even if rj = 
so in the future, we need not insert provisos. 

For use in the remainder of this proof, we shall divide T d L into disjoint cubes 
Ci, . . . Cj, . . . (half open/closed etc.) of diameter a. Since we are only pursuing 
the limit of 9[ L \ it may just as well be assumed that 2a divides L. We use the 
notation 

W/Wlhcj) ■■= [ \fW\da: 

JCj 

and similarly for other local norms. 

Our first substantive claim is as follows: 
Let e denote a small number which is of order unity independent of L (the nature of 
which is not so important and will, to some extent, be clarified below) and suppose 
that for all j, 

\\v\\LH Cj )<^ (68) 

Then, for all 9 is sufficiently small, for all L under discussion, we have that r\ = 0. 
To see this we let x S Cj be in the support of rj + so, ostensibly, we have 

l + i]^(x)=e e l . (69) 

However, due to the finite range of V, the integration actually takes place on only 
the cubes in the immediate vicinity of Cj so that 

1 + V + (x) < exp[6V max Y, IMIli(CV)] < e eK,,axDl (70) 

where f ~ j means that Cf D Cj ^ and D\ = D\{d) is the number of j' such 
that j' ~ j. This implies an L°°-bound on r] + which is (a small number) of order 
unity. We run a similar argument for ry~ - only now we have to contend with k: 

1 _ v - > e -« e -(V«fli ( 71 ) 

i.e., 

yf < K + eVmuDL (72) 

Thus we have an L°° bound on the full r] which implies, at this stage - since h and 
are supposed to be small - an improved bound on ||j7||ii(c^) i n ai l cubes Cj. Let 
us continue the process. Suppose that for all j 



MIli(co < ( 73 ) 
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where <fi = (p(9,h) represents the latest improvement. Then 

1 + Halloo <exp[Di\W0] (74) 

- so that | |r? + 1 |oo < DiV max 4>9 - and 

IMU < DiV m ^9 + K (75) 

which, at least for a while, represents an improvement on the various L 1 -norms. 

The procedure is no longer beneficial if cj) is on the order of k. E.g. we may 
stop when k > 4>[[2a]~ d — 9V max Di] - where it is assumed that 9 is small enough 
so that the coefficient of <j> is positive. We arrive at an overall bound: 

IMloo < C a K < c b 9h 

where the c's are constants of order unity independent of L and 9 (provided that the 
latter is sufficiently small). Clearly, for 9 small enough, this cannot be consistent 
with ||??po||i = h unless h = 0. 

Thus, we are done with the proof unless there are bad blocks where the local 
L 1 -norm of rj is in excess of e9~ 1 . In fact, we will present an additional hierarchy 
of bad blocks. The above blocks will be the core blocks which will be denoted by 
C. We shall define blocks B n , n = 0, 1, . . . s by 

B n = {c j | e r<|M| L1(c . ) <r} (76) 

and (unfortunately) other bad blocks 

B' n = {C j \ee n - 1 >\\ V \\ LHc . ) >e n }. (77) 

Thus it is seen that our e should be small enough to absorb various constants which 
crop up but large compared to the assumed value of 9. The hierarchy of these 
sets stops when the L 1 -norm is comparable to k - pretty much as in the previous 
argument. Here we shall say that s is defined so that blocks outside the hierarchy 
have local L 1 -norm of rj less than a constant Qq times k with Qq to be described 
shortly. Such blocks will be informally referred to as background blocks. 
We order the hierarchy in the obvious fashion: 

... B' n yB n yB' n+1 y ... (78) 
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with every element of the hierarchy considered to be above the background blocks 
and below the core blocks. 

Our next claim is that for any block in the hierarchy there must be a neighboring 
block which is further up the hierarchy. Indeed suppose not and that, e.g., Cj £ B n . 
If all neighbors - that is to say all blocks Cy with j' ~ j - were only at the level 
B n or below, we would have, for x £ Cj 

l + rj + (x) < e eDlKm " e " (79) 

and 

rj-(x)<9D 1 V max 9 n + K. (80) 

This, for most cases, implies that | \r]\ \l°°(Cj) * s °f trie order 9 n+1 which precludes 
Cj £ B n ; a similar derivation applies to the B^. At the very bottom of the hier- 
archy, the same situation holds with any reasonable choice of Qo which is of the 
order of unity. 

The implication of the preceding claim is that each block in the hierarchy is 
connected to the core by a path whose length does not exceed the order of its the 
hierarchal index. 

Our next claim is that (under the hypothesis of non-triviality) the vast majority 
of the L -norm of rj is carried by the core and its immediate vicinity. First, let us 
estimate |C|, the volume of the core. We argue that 

\C\-e-\<hL d (81) 


since the right side is the full L 1 -norm of rj and the left side represents the minimal 
contribution to this effort on the part of the core. Thus the core volume fraction is 
the order of 9h which, we remind the reader is supposed small. 

Now, by the connectivity property of the hierarchy, Eq.(81) can be translated 
into an estimate on the volume of the sets B„, ~B' n . Indeed, we may write 

|B n | < \C\D 2 n d (82) 

where D 2 = D 2 (d) is another constant. 

The above two estimates are sufficient to vindicate the claim at the beginning 
of the paragraph containing Eq.(81). We denote by C* = C U Bo U B' which 
we call the extended core. Turning attention to the complimentary set we have, for 
B n : 

\rj\dx < 9 n \C\D 2 n d [2a]- d < h9L d -^ ■ 9 n n d D 2 [2a]- d (83) 
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Summing this over all n starting from n = 1 we get a contribution (for all 9 suffi- 
ciently small) which is bounded by Ke~ l 6 2 hL d for some constant K. Similarly 

/ \ V \dx < e8 n - 1 \C\D 2 n d [2a]- d (84) 

whence the total contribution from the primed portion of the hierarchy on the com- 
pliment of the extended core to is no more than K'9hL d for some constant K' 
Note that these are small compared to the purported total of hL d . Finally, the con- 
tribution from all the background blocks is surely no more than Qo[2a]~ d k[L / a] d 
and we recall that k itself is bounded above by the order of 9h. So, in summary, 
we arrive at 

/ rj + dx > (1 - g9)hL d (85) 

for some constant g which is independent of L and 9 (for 9 sufficiently small). 

With the preceding constraint in mind, let us bound from below the relative 
entropy 

S(p)-S(po)= I po(l + 7? + )log(l + r? + )+ / p (l-??-)log(l -??~)- (86) 

The second term may be bounded below by — \h. As for the former, since the 
function is always positive, we may restrict attention to the set C*. As is not 
hard to show, the contribution from C* is larger than that of the function which is 
uniform on C* and has the same total mass. As a result: 

S(p) - S(po) > -\h + S (1 + -Ul - 9 9)hL d ) log(l + y^(l " 99)hL d ). 

(87) 

As is not hard to see, (and is intuitively clear) this is decreasing in | C* | - meaning 
we may substitute the upper bound based on Eq.(81): |C*| < L d G9h with G 
(oc e -1 ) another constant of order unity independent of L and 9. All in all, 

S{p) - S( P0 ) > -h + G8h(l + { ^r^\) log(l + ^~)- (88) 

By contrast we have, from Eq.(63), that S{p) — S(po) is less than a constant 
times Oh. This along with Eq.(88) implies that h = or, assuming that ft / Oa 
strict lower bound on 9. Either of these conclusions allows us to infer the desired 
result. □ 
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3.3 Catastrophic Behavior 

We conclude with some examples of what can go "wrong" if the criterion of ther- 
modynamic stability is violated. For the benefit of these final results we shall vio- 
late condition-K by assuming the existence of a compactly supported probability 
density Pq(x) such that 

V(x - y)p Q {x)p & (y)dxdy = -u (89) 
with uo > 0. First, some specific instances: 

Proposition 3.3. IfV satisfies either of the following then it violates condition— K 

(a) The interaction V satisfies f Rd V(x)dx = — Vq < 

(b) For some A < 1, in a Aao neighborhood of the origin, V integrates to +cq 
while for Aao < \ x \ < 2ao, V(x) is bounded above by —vq where 

«o(l - X d ) > cq. 

Proof. In the first case, we choose 



!,©0) = X\x\<ta-T7-u 



1 

1 tN' 



where 7 is a geometric constant. It is noted that 

Km g{at) d / V(x - y)p^ Q (x)p itQ (y)dxdy = -v (90) 

so the result follows for £ sufficiently lai - ge. 

As for the second case, we simply use p&(x) = ^^ d X\x\<a - m performing 
the integration 

/ V(x - y)p & (x)p Q (y)dxdy 

JR d xR d 

and ignoring the positive contribution from \x — y\ < Aao the result would be not 
more than — vq. For each x we must cut out a ball of radius (no more than) Aao 
around each point of the integration and insert a corresponding factor of (no more 
than) Co- The result is no more than — («o(l — A rf ) — cq) and we are done. It is 
noted that this latter result applies immediately to the case where V is negative in 
a deleted neighborhood of the origin. □ 
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Theorem 3.4. For a potentials that violate condition— K via Eq.(89), the V-McK 
system exhibits non— physical scaling. In particular, for L sufficiently large, 

9 T (L) < rL~ d 

for some r > 0. 

Proof. Recalling Eq.(49) we have for any L 

By contrast, if we abide by the recommended density we obtain: 

where it is noted that S(p & ) < oo by hypothesis and by the restrictive nature of the 
support of p & , it is independent of L. But then, as soon as — ^9uoL d + S(pq) < 
log L d + ^6go it must be that 6 < 9 T . This implies the stated bound. □ 
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